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ABSTRACT 

We used the 1.4 GHz NVSS to study radio sources in two color-selected QSO samples: a volume- 
limited sample of 1313 QSOs defined by M\ < —23 in the redshift range 0.2 < z < 0.45 and a 
magnitude-limited sample of 2471 QSOs with m r < 18.5 and 1.8 < z < 2.5. About 10% were detected 
above the 2.4 mJy NVSS catalog limit and are powered primarily by AGNs. The space density of 
the low-redshift QSOs evolves as p oc (1 + z) 6 . In both redshift ranges the flux-density distributions 
and luminosity functions of QSOs stronger than 2.4 mJy are power laws, with no features to suggest 
more than one kind of radio source. Extrapolating the power laws to lower luminosities predicts the 
remaining QSOs should be extremely radio quiet, but they are not. Most were detected statistically 
on the NVSS images with median peak flux densities S'pfmJy beam -1 ) pa 0.3 and 0.05 in the low- 
and high-redshift samples, corresponding to spectral luminosities log[Lx.4GHz(W Hz -1 )] « 22.7 and 
24.1, respectively. We suggest that the faint radio sources are powered by star formation at rates 
M ~ 20M© yr -1 in the moderate luminosity (median (Mi) « —23.4) low-redshift QSOs and M ~ 
500M Q yr -1 in the very luminous ((-Mi) ~ —27.5) high-redshift QSOs. Such luminous starbursts 
[(log(LjR/L©)) ~ 11.2 and 12.6, respectively] are consistent with "quasar mode" accretion in which 
cold gas flows fuel both AGN and starburst. 

Subject headings: galaxies: active — galaxies: starburst — quasars: general — radio continuum: galaxies 



1. INTRODUCTION 

The radio emission from quasi-stellar objects (QSOs) 
has puzzled astronomers for nearly 50 years. The first 
QSOs to be recognized contain very strong radio sources 
similar to those powered by the active galactic nuclei 
(AGNs) of the most luminous radio galaxies, but ~ 90% 
of bright optically selected QSOs are fainter than 1 mJy 
at 1.4 GHz. Also, the flux-density distribution of radio- 
loud QSOs is quite different from that of radio galax- 
ies, and might even be bimodal. These differences led 
to longstanding debates about whether QSOs experience 
the same cosmological evolution as radio galaxies and 
whether there are two distinct types of radio sources in 
QSOs. In this paper, we use direct detections of radio 
sources stronger than 2.4 mJy and statistical detections 
of fainter sources in two large samples of QSOs to argue 
that (1) radio-loud QSOs are similar to powerful radio 
galaxies and do evolve, despite their differing flux-density 
distributions, and (2) radio-quiet QSOs really are differ- 
ent from radio-loud QSOs because most of their radio 
emission is powered by starbursts instead of AGNs. Since 
the problem of radio emission from QSOs has such a long 
and complicated history, we begin with a short review to 



set the stage and define common terms. 

Quasars, QSGs, and QSOs 
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(|i"963h . and 
the first quasi- 
stellar radio sources (quasars or QSSs) by optically 
identifying redshifted "stars" with radio sources whose 
luminosities are "not markedly different from other 
known strong radio sour c es lik e Cygnus A or 3C 295" 
(jGreenstein fc Matthews! 119631 ) in elliptical galaxies. 
The first o ptically s e lected quasi-stellar galaxies (QSGs) 
found by ISandagd (|1965| ) "resemble the quasi-stellar 
radio sources in many optical properties, but they are 
radio quiet" in the sense that only about one QSG 
per thousand brighter than m pg = 19 appeared in the 
relatively insensitive all-sky radio catalogs of the time 
(S > 9 Jy at 178 MHz). All of the first quasars and 
QSGs appear starlike because their AGNs are far more 
luminous than their inferred host galaxies, assuming 
their redshifts are cosmological. The deliberately 
vague name "quasi-stellar object" (QSO) still used 
today for both qu a sars a nd QSGs was introduced by 
iHoyle &: Burbidgel (|1966| ) because they doubted the 
redshifts of the QSGs were cosmological and so "at this 
stage we do not think enough is known about their 
nature for a definitive name to be chosen." 

Although QSGs were called "a m ajor new constituent 
of the universe" by Sandagej (|1965l ). they are not a dis- 
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tinct constituent, only the luminous tail of a continuous 
distribution that includes Seyfert gal axies whose AGNs 
do no t overwhelm their starlight. Sc hmidt fc Greenl 
(1983) formally distinguished them from Seyfert galaxies 
by requiring QSOs contain AGNs brighter than Mb = 
—23, the absolute magnitude of the brightest galaxies 
of stars calculated using Hq = 50 km s _1 Mpc -1 and 
g<o = 0.1. In the modern flat ACDM cosmology with 
H = 71 km s _1 Mpc -1 and D, m = 0.27, the bright- 
est galaxies have Mb ~ —22.2, so the spirit of the 
ISchmidt fc Greenl (| 19831 ) criterion (QSOs contain AGNs 
brighter than any host galaxy of stars) would change the 
QSO cutoff to Mb ~ —22.2. However, like most authors, 
we fol low the letter (M B = -23) of the ISchmidt fc Greenl 
(1983) criterion. Thus our least-luminous QSOs are ac- 
tually about twice as luminous as the brightest galaxies. 
The value of this cutoff matters because such luminous 
QSO s are rarely found in spiral galaxies (jDunlop et al.l 
120031) . 

1.2. QSOs Versus Radio Galaxies 

The cumulative radio flux-density distribution of the 
strongest radio galaxies rises slightly faster than the Eu- 
clidean rate 

/>oo 

JV(> S) = / n{s)als oc S~ 15 , (1) 
Js 

where A^(> S) is the cumulative number of sources per 
steradian stronger than S and n(S)dS is the differential 
number per steradian with flux densities S to S + dS. 
If the same were true of bright QSOs and there is one 
radio-loud QSO stronger than 10 Jy among 1000 opti- 
cally selected QSOs, one might hope to find 30 radio 
sources stronger than 1 Jy and hundreds above 0.1 Jy. 
From this perspective, the first directed searches for ra- 
dio emission from o ptically selected QSOs were surpris - 
ingly unsuccessful (jKellermann fc Paulinv-Tothl [1966). 
Radio sources were detected in only about 10% of 
bright QS Os even after decad e s of increasingly sensitive 
searches dKatgert et all 119731: IStrittmatter et al.l 119801: 
IMiller et al.l "fl990 N ) for sources as faint as 0.01 Jy at 
v = 1.4 GHz, orlog[L l/ (WHz -1 )] w 24 at z « 0.2. These 
low detection rates imply that the flux-density distribu- 
tion of optically selected QSOs rises much more slowly 
and may even be bimodal. They might also indicate 
that QSOs and radio galaxies do n ot experience the sa me 
amount of cosmological evolution (Mi ller et alJ lT990). 

1.3. Are There Two Distinct Populations of QSOs? 

IStrittmatter et all (|1980l ) were the first to report a bi- 
modal QSO flux-density distribution, with an observed 
peak near S = 1 Jy and an inf erred peak at S < 1 mJy, 
albeit in a very small sample. IMiller et al.l (|1990f ) noted 
a bimodal luminosity distribution implied by the bi- 
modal flux-density distribution of their larger sample 
of QSOs in the redshift range 1.8 < z < 2.5. They 
found nine radio sources with v = 5 GHz spectral lu- 
minosities log[-L„(WHz -1 )] > 26, 96 upper limits be- 
low logfL^WHz -1 )] rs 25, and nothing in between. 
The term "bimodal" is perhaps too strong: a bimodal 
distributi on has two p eaks. Neither IStrittmatter et all 
(|1980D or IMiller et all (|1990l) actually detected the low- 
luminosity peak; its existence was only inferred from the 



large number of nondetections whose radio flux densities 
or luminosities would have to be spread smoothly over 
many orders of magnitude to not yield a second peak. 
For a recent analysis of b imodality in light of m odern 
optical and radio data, see iBalokovic etafl (|2012T >. 

A bimodal flux-density or luminosity distribution sug- 
gests there may be two distinc t populations of r adio 
sources in QSOs. For example, IMiller et al.l (|1990l) as- 
sociated their high-luminosity peak with AGNs in ellip- 
tical galaxies similar to the hosts of radio galaxies and 
the inferred low-luminosity peak with weaker AGNs in 
spiral galaxies. However, any nonuniformity or feature 
in these distributions is sufficient to suggest two distinct 
populations. Bimodality is sufficient but not necessary, 
so we suggest that debates about whether a clearly vis- 
ible feature in a distribution makes it truly bimodal are 
not relevant to the fundamental question "Are there two 
distinct populations of QSOs?" 

1.4. Further Complications 

Three factors have complicated investigations of QSO 
radio luminosity functions: 

(1) The radio/optical luminosity ratio of QSOs might 
be a better measure of radio loudnes s than absolute ra- 
dio luminosity alone. ISchmidt] ()1970D proposed that the 
radio luminosity function "depends on the optical lumi- 
nosity in such a manner that the ratio has a universal 
distribution function" of the form 

$(F ,F I ) = <Z>(F )*(F T /F ) . (2) 

Equation[2]implies a linear correlation between radio and 
optical luminosities, so the radio/optical ratio 

R = F T /F (3) 

may b e a better indicator of radio activity. IPeacock et al.l 
(1986) countered that almost no radio-loud QSOs are op- 
tically fainter than M B = -24 (for H = 50 km s Mpc -1 
and go = 0.1, equivalent to Mb > —23.2 for to- 
day's ACDM model), so that an apparent correlation 
in the data may only reflect differences between two 
populations of QSO host galaxies, optically luminous 
ellipticals and optically fainter spirals. Using FIRST 
(jBecke r. White, & Helfand 1995) direct a nd statistical 
detect ions of QSOs in a very large sample, iWhite et al.l 
(|2007l ) found a clear but nonlinear correlation of radio 
and optical luminosities F r oc F° s5 . However, Appendix 
C of llvezic et ail (2002) argues that selection effects are 
responsible for such apparent correlations. 

(2) Flux-density boosting in rclativistic [y 2 = (1 — 
v 2 /c 2 )~ 1 3> 1] radio jets may significantly broaden both 
the radio flux-density and the R distributions of QSOs, 
whose op tical emission does not app ear to be strongly 
beamed (Sch euer fc Readheadl Il979). Thus radio jets 
in optically selected QSOs should be randomly oriented 
relative to our line-of-sight, but only the small frac- 
tion m 7~ 2 of QSOs whose jets are aligned near our 
line-of-sight become radio-loud quasars. The geometry 
of relativistic beaming tends to flatten the differential 
count of QSO radio sources from n(S) oc S~ 5 ^ 2 towards 
n(S) oc S-^ 3 [or N(> S) oc S' 1 ^ 3 } for intrinsically 
steep-spectrum (a = d\nS/dhiv ~ —1) sources. 

(3) Like the optical emission from a QSO, the ra- 
dio emission is the sum of two known components, 
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one powered by the AGN and the other by the host 
galaxy of stars. Massive stars born within the past 
~ 10 8 years and relativistic electrons accelerated by 
their supernova remnants power the radio emission 
from the host galaxy, so a "red and dead" ellipti- 
cal galaxy with no recent star forma tion can be ex- 
tremely radio-quiet ([Walsh et al.lll989|) . A number of 
features can be used to distinguish between these two 
types of radio source. The most powerful nearby star- 
burst is in the ultraluminous infrared galaxy Arp 220, 
whose spectral luminosity is log[Li.4GHz 

(WHz" 1 )] w 

23.5. Although the western nucleus of Arp 220 ap- 
pears to contain an obscur ed AGN with X-ray lumi- 
nosity L x ~ lO^ergs" 1 ([Rangwala et al.ll201lD . the 
compact radio component in the western nucleus ac- 
counts for less than 25% of the integrated 1 . 4 GH z 
radio flux density (|Mundell. Ferruit. fc Pedlar! 1200 lh . 
so the Arp 220 starburst luminosity is not less than 
log[Li.4GHz(W Hz -1 )] w 23.4. A starburst-powered ra- 
dio source is roughly coextensive with the starburst and 
is not brighter than T b ~ 10 5 K at 1.4 GHz, has the 
steep spectrum a ~ —0.7 at centimeter wavelengths typ- 
ical of optically thin synchrotron r adiation, obey s the 
far-infrared (FIR) /radio correlation (|Condonlll992l) . and 
is not strongly variable. In contrast, an AGN can pro- 
duce a radio source several orders of magnitude more 
luminous than log[Li.4GHz (WHz -1 )] ~ 24, its radio 
jets and lobes can extend far outside the host galaxy 
of stars, and its radio core can be s mall enough to 
vary significantly on time scales of years (Barvainis et al.l 
2005). The radio spectrum of an AGN can be nearly 
flat (a ~ 0) over a wide frequency range, either from 
a synchrotron self -absorption "conspi racy" among com- 
pact components ([Cotton et al.lll980T ) or from optically 
thin free-free radiation by a hot T > 1 7 AGN disk wind 
' Blundell fc Kuncid I2007T ) . However , iSteenbrugge et al.l 



20111 ) found that "for 20 out of 22 PG quasars stud- 



ied, free-free emission from a wind emanating from the 
accretion disc cannot mutually explain the observed ra- 
dio and X-ray luminosity," eliminating free-free emission 
as the cause of the radio-quiet peak in the luminosity 
distribution of QSOs. 

Nearly all low-redshift (z < 0.5) radio galaxies 
seem to be "red and dead" ellipticals whose radio 
sources are powered primarily by their AGNs. By 
analogy, low-redshift QSOs whose radio sources are 

(1) more luminous than log[Li.4GHz (WHz -1 )] - 24, 

(2) have lower FIR/radio flux-density ratios than star- 
burst galaxies, ( 3) are either larger than ~ 20 kpc 



or ver y compact dFalcke et al1ll996t iBlundell fc Beasleyl 
119981: lUlvestad et all I2005D. (4) are strongly variable 



( Barvainis et al.ll2005 ) or (5) have flat radio spectra prob- 
ably reside in elliptical galaxies. This analogy is sup- 
ported by hig h-resolution HST images of nearby QSO 
host galaxies ([Bovce et all 119981: iDunlop et all ( 2003). 
QSOs close to the original iSchmidt fc Greenl ( 19831 ) lu- 
minosity limit often obey the FIR/radio correlation, 
may contain s t ar burst-dominated radio sources (e.g. , 
Barvainisl [19901: ISopp fc Alexanderl fl99lt IBarthell [2001 
de Vries et al.l I2007D. and may be in dis c or merging 
galaxi es ([Bovce et al.l [19 98) . In contrast. IDunlop et al.1 
(2003) found that QSOs having nuclear luminosities V < 
-23.5 (assuming H = 50 km s" 1 Mpc" 1 , or V < -22.7 



in modern ACDM cosmology) are virtually all massive el- 
lipticals with no large-scale star formation. Weak radio 
emission correlated with the strength of the 4000 A break 
indicates that star formation is qui te common in the hos t 
galaxies of low- luminosity AGNs ([de Vries et al.l [2007). 
We conclude that the radio properties of QSO samples 
are likely to depend on the absolute magnitude cutoff, 
and samples containing low-luminosity AGNs should not 
be used to infer the radio properties of true QSOs. 

The 1.4 GHz flux density of an ultraluminous star- 
burst galaxy like Arp 220 would be only S « 0.4 mJy 
at z « 0.5 and S ~ 0.015 mJy at z 2, so excep- 
tionally sensitive observations are needed to character- 
ize the starbur st radio emission from individual QSO 
host galaxies. iKellermann et all (|1989f ) detected most 
of the low-redshift (z < 0.5) Palomar Bright Quasar 
Survey (PBQS) QSOs at 5 GHz and found a bimodal 
spectral luminosity distribution (not function) peaking 
at log[L 5 (WHz -1 )] ~ 26 and ~ 22. The low-luminosity 
peak for QSOs spans the range 0.1 < R < 1 with a me- 
dian (R) ~ 0.3. iMiller et al.l |l993l ) used high-resolution 
radio images and [O III] spectra of the PBQS QSOs 
as evidence for a circumnuclear starburst contribution 
to the low-luminosity peak. Recent results from signifi- 
cantly more sensitive Karl G. Jansky Very Large Array 
(VLA) observations of 179 SDSS (Sloan Digital Sky Sur- 
vey) QSOs in the reds hift range 0.2 < z < 0.3 found 
a low-luminosity peak (|Kimball et al.|[2~011[ ) centered on 
log[L 6 (W Hz -1 ) « 22.4 and R « 0.3 at v = 6 GHz. 
The radio sources contributing to this peak have a me- 
dian spectral index a = —0.7 consistent with starburst 
radio emission. Unfortunately their median flux den- 
sities are only (5) ~ 0.13 mJy at 6 GHz and (S) fts 
0.33 mJy at 1.4 GHz so making foll o w-up observations 
will be difficult. The iKimball et "all (|2011l ) radio lumi- 
nosity function suggests that most radio sources with 
log[L 6 (W Hz" 1 )] > 23 at 6 GHz (or log[Li. 4 (W Hz" 1 ) > 
23.4 at 1.4 GHz) in low-redshift QSOs are powered pri- 
marily by AGNs, no t starbursts. 

IKellermann et all (|1994f ) reobserved "radio quiet" 
QSOs and low-luminosit y AGNs from their BQS sample 
(IKellermann et alJ lT989) at 5 GHz with 18" and 0'.'5 res- 
olution and resolved a number of compact cores and ex- 
tended lobes similar to those in radio-loud quasars, sug- 
gesting that central engines exist in galaxies and quasars 
over a wide range of radio luminosity. They noted, how- 
ever, that their resolved "radio quiet" QSOs had radio 
luminosities well in excess of that expected from the un- 
derlying galaxy. 

Several high-resolution radio observations have demon- 
strated that most mJy-level radio sources in QSOs con- 
tain radio components so compact that they must be 
AGN-powered. Three flat-spectrum (a > —0.5) "ra- 
dio intermediate" (defined by 10 < R < 250 at v = 5 
GHz) QSOs were found to have brightness temperatures 
T b > 10 9 K ([Falcke et alJ[l99l . much too high for 
starburst radio emission. However, only one of them 
(PG 1309+355 at z = 0.184) satisfies our M B < -23 
criterion for a QSO, and its 5 GHz radio luminos- 
ity log[L 5 (W Hz" 1 )] rj 24.6 is an order of magnitude 
higher than the m ost luminous starbu r st yet found in a 
low-redshift QSO. IBlundell fc Beaslevl (fl9M ) imaged 12 
"radio-quiet quasars" with the VLBA and found strong 
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evidence for jet-producing central engines in eight of 
them. However, four of the "quasars" are low-luminosity 
AGNs with M B > -23 (Mrk 335, III Zw 2, Mrk 705, 
and II Zw 171). Five of the eight true QSOs were 
detected by the VLBA, but all five have radio lumi- 
nosities in the range 23.8 < log[Li. 4 (W Hz -1 ) < 26.2] 
high eno ugh that we expect a significant AGN contri- 
bution. lUlvestad et al.l ([20051 ) considered a sample of 
11 "radio-quiet quasars" having total 5 GHz flux den- 
sities of 4 mjy or greater. Three have no compact 
cores stronger than 1 mjy and were discarded from the 
sample. New VLBA images were made of five, one of 
which (UGC 09412) is a lo w-luminosity AGN. Two o f 
the true quasars are in the iBlundell fc Beaslevl (fl998h 
sample above. The remaining QSOs (J0804+6459 at 
z = 0.148 and UM 275 at z = 2.137) have radio luminosi- 
ties log[L 5 (W Hz -1 )] = 23.8 and 26.0, respectively, and 5 
GHz radio/optical flux ratios R — 3 and 22, respectively, 
much higher than we expect for radio sources powered 
only by starbursts. We believe their conclusion "Thus, 
the most likely explanation for their radio properties is 
simply the radio-quiet quasars are similar to their 'tra- 
ditiona l' radio-loud cousins but have less powerful radio 
jets." ([Ulvestad et al.l [2005]) applies only to QSOs con- 
taining radio sources stronger than starbursts produce. 
We conclude that, unfortunately, VLB observations are 
simply not sensitive enough (a <C 0.1 mjy) to address 
the question of starburst radio emission in large samples 
of radio-q uiet quasars. 

Finally, Barvaini s et al.l ([2005) compared radio vari- 
ability in samples of radio-quiet (defined by having R < 
3), radio- intermediate, and radio- loud quasars. They 
found that radio variability is independent of radio loud- 
ness. Of the 11 "radio-quiet quasars" monitored, four are 
low-luminosity AGNs (Mrk 335, Mrk 1148, UGC 11763, 
and PGC 070504). Four of the seven true QSOs are vari- 
able radio sources and three are not. Among the variable 
QSOs, only PG 0052+251 at z = 0.155 is sufficiently ra- 
dio quiet (log[L 8 (W Hz -1 )] = 22.6 and R ~ 0.3) that a 
starburst might have contributed significantly to its ra- 
dio emission. This single source is consistent with the 
6 GHz radio lumino sity functions shown in Figure 6 of 
iKimball et al.l ([201 ll) . which suggests that about one in 
four QSO radio sources with log[L 6 (W Hz -1 )] = 22.6 
should be AGN-dominated. We believe their conclu- 
sion "...the radio emission from radio-quiet quasars orig- 
inates in a compact structure intimately associated with 
the active nucleus. The alternative hypothesis, that the 
emission from radio-weak quasars is from a starburst, 
is ruled out." ([Barvainis et al.ll2005h has been justified 
only for radio sources stronger than log[Ls(W Hz)] ~ 24 
and R > 1 in QSOs, not for the weaker radio sources 
expected from starbursts. 

1.5. Outline of This Paper 

QSOs have such broad radio luminosity functions that 
their radio emission can be characterized only via large 
unbiased samples of radio-detected QSOs. Sensitive 
directed searches for radio emission from QSOs have 
been impeded by the small numbers of radio detec- 
tions. A typical VLA observing program might target 
N ~ 100 QSOs. Detecting statistically useful numbers 
N ~ 100 ^ V 1 / 2 of radio sources requires either a de- 



tection rate near 100% or a much larger QSO sample if 
the detection rate is low (e.g., N ~ 1000 if the detection 
rate is only ~ 10%). 

Furthermore, nearby QSOs in which starburst-powered 
radio sources can be detected individually have such a 
low sky surface density that a statistically useful sam- 
ple can be obtained only from an optical survey cover- 
ing more than a steradian with uniformly accurate mul- 
ticolor photo metry. The SPS S DR7 (Data Release 7) 
QSO sample ([Schneider et aLllMToh covers f2 ~ 2.66 sr 
at Galactic latitudes \b\ > 30° and satisfies that require- 
ment. T his paper is based o n NRAO VLA Sky Survey 
(NVSS) (|Condon et al.l fl998h 1.4 GHz flux densities of 
(1) a volume-limited sample of 1313 DR7 QSOs with 
0.2 < z < 0.45 (Sec. [2]) and (2) a magnitude-limited sam- 
ple of 2471 DR7 QSOs in the redshift range 1.8 < z < 2.5 
(Sec, E} . 

Above the NVSS catalog limit S « 2.4 mjy beam -1 , 
the detection rate in each sample is only about 10%, 
but the QSO samples are large enough to yield hundreds 
of NVSS detections. In addition, most of the remain- 
ing QSOs can be detected statistically on the NVSS im- 
ages at levels S ~ 0.1 mjy beam -1 . The main result of 
the direct detections is that the differential luminosity 
function p m (L) of strong radio sources powered by the 
AGNs in QSOs is an extremely flat power law, so flat 
that extrapolating to lower luminosities predicts essen- 
tially no statistical detections at the S ~ 0.1 mjy level. 
Our high statistical detection rate implies a "bump" 
in the radio luminosity function of QSOs in the range 
log[L(WHz -1 )] - 23 to 24, depending on redshift. We 
attribute this bump to radio emission powered by lumi- 
nous starbursts in the host galaxies of most bright QSOs 
(Sec. P. 

2. THE VOLUME-LIMITED 0.2 < Z < 0.45 QSO SAMPLE 

The SDSS DR7 QSO catalog (jSchneider et al.ll27)Toh is 
complete to i = 19.1 mag over a solid angle = 2.66 sr 
around the North Galactic Pole. It contains the small 
sample of 179 color-selected QSOs defined by Mi < —23 
in the narrow redshift range 0.2 < z < 0.3 studied 
by IKimball et al.l ([201 ll ) and the larger sample of 1313 
QSOs in the wider redshift range 0.2 < z < 0.45 dis- 
cussed here. Note that these magnitudes were calculated 
for an Ho = 71 km s Mpc -1 ACDM cosmology, so our 
QSO cut off corresponds to Mb ~ —23.8 for compari- 
son with iSchmidt fc Greenl ([1983D and should discrimi- 
nate against the low-luminosity QSO s that can appear 
in spiral galaxies ([Dunlop et al.ll2003l ). Both samples are 
strictly volume limited because any Mi < —23 QSO with 
z < 0.45 is brighter than i = 19 mag. There are no 
more color-selected QSOs with 0.2 < z < 0.45 to be 
found in that 21% of the sky, so no observation made 
at this time and place in the universe can reduce our 
statistical errors by more than about a factor of two — 
our sample variance is approaching the "cosmic variance" 
limit set by the finite observa ble volume of the universe 
([Kamionkowski fc Loeblll997l ). 

2.1. Sources in the NVSS Catalog 

The entire DR7 area is covered by the NVSS, whose 
source catalog is complete for statistical purposes above 
a peak flux density S p ~ 2.4 mjy beam -1 at 1.4 GHz. 
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In the redshift range 0.2 < z < 0.45 the 45 arc- 
sec FWHM (full width between half-maximum points) 
beam of the NVSS spans 150 to 250 kpc, more than 
enough to cover the host galaxy of stars and most of 
the extended radio emission powered by an AGN. Con- 
sequently the NVSS sees most QSOs as point sources 
and the NVSS integrated flux-density limit is S « S p = 
2.4 mjy. For radio sources with spectral index a ~ —0.7, 
the corresponding v = 1.4 GHz (in the source frame) 
NVSS spectral-luminosity detection limits range from 
log[L !/ (WHz -1 )] = 23.4 at z = 0.2 to 24.2 at z = 0.45. 
Thus most of these sources are too luminous to be pow- 
ered primarily by starbursts. 

To ensure completeness even for radio sources larger 
than the NVSS beam, we examined all NVSS catalog 
entries within 120 arcsec of each QSO. We rejected com- 
ponents centered more than a few arcsec from the op- 
tical position unless their morphologies in the NVSS 
and/or higher-resolution (5.4 a rcsec FWHM) FIRST 
(|Becker. White, fc Helfaiidl [19951) images plausibly link 
them with the QSO. Unresolved or symmetric NVSS 
sources are easily identified with QSOs by radio-optical 
position coincidence alone, but asymmetric multicom- 
ponent sources can be difficult. For example, the left 
panel of Figure [1] shows two marginally resolved NVSS 
components whose centers straddle the optical position 
(cross) of SDSS J091205. 16+543141.2. Two lobes and 
a core are clearly resolved in the FIRST image (right 
panel), confirming the otherwise uncertain NVSS iden- 
tification. Because the NVSS has low resolution, it has 
high surface-brightness sensitivity. The two NVSS com- 
ponents contain essentially all of the QSO flux density, 
23.8 mjy, while the three FIRST components add up to 
only 6.6 mjy. This example illustrates how the NVSS 
and FIRST complement each other for identifying ex- 
tended radio sources and measuring their flux densities. 

We also found that the FIRST peak flux densities of 
QSOs fainter than S « 10 mjy are usually close to the 
NVSS flux densities, and very few faint NVSS sources are 
not in the FIRST catalog. This indicates that most of the 
fainter radio sources produced by these QSOs have half- 
power linear sizes < 5 kpc and the likelihood that the 
NVSS has resolved out any faint sources is small. Unfor- 
tunately, radio sources powered by starbursts and and 
by low-luminosity jets from AGNs are both consistent 
with this size limit; sub-arc second resolution is needed to 
distinguish between them ([Condon et al.lll991ft . FIRST 
probably does not resolve out many sources in the 1 < 
S'(mJy) < 2.4 range either, but we have not included 
the FIRST flux densities for these sources in our sta- 
tistical analysis because they are systemat ically low by 
abou t 0.25 mjy owing to "snapshot" bias (I White et"aLl 
120071) . Also the FIRST source count (jWhite et alj|1997ft 
indicates a rapid fall in differential completeness below 
S~2 mjy. Although the NVSS image noise is higher, 
the snapshot bias is lower, so we chose to study all sources 
fainter 2.4 mjy statistically via their amplitude distri- 
bution on NVSS images (see Section [2751 for a complete 
discussion of snapshot bias in FIRST and NVSS) instead. 

A number of stronger sources are clearly resolved by 
the NVSS, but most have brightnesses well above the 
2.4 mjy beanT 1 NVSS catalog limit. Thus the NVSS 
integrated flux densities of radio-loud QSOs should be 
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Figure 1. NVSS (left panel) and FIRST (right panel) images 
of SDSS J091205. 16+543141.2 (cross). Contours are separated by 
factors of ±2 1//2 starting at 1 mjy beam -1 (NVSS) and 0.5 mjy 
beam" 1 (FIRST). 

accurate and the NVSS should not have missed any 
radio-loud QSOs. There are 163 NVSS detections with 
0.2 < z < 0.45 (Table [J), of which only 37 are in the 
narrow redshift range 0.2 < z < 0.3 that was also cov- 
ered by our more sensitive d irected VLA observations at 
6 GHz ([Kimball et alJl20TTh . 

2.2. The 1.4 GHz Flux-density Distribution of QSOs 
Stronger Than 2.4 mjy 

Let n(S)dS be the differential number of sources per 
steradian in the flux-density range S to S + dS. Then 
the brightness-weighted count S 2 n(S) is proportional to 
the contribution dT^ of sources in each logarithmic flux- 
density range to the sky background temperature: 



dT h 




'ln(10)c 2 " 


[d\og{S)_ 




2k B v 2 



S 2 n(S) . (4) 



The brightness-weighted count of all extragalactic 
sources at 1.4 GHz (mostly radio galaxies) is shown by 
the upper curve in Figure [3J The lower curve in Fig- 
ure [2] shows that (1) the 1.4 GHz flux-density distribu- 
tion of the 0.2 < z < 0.45 QSOs is a smooth power law 
above 2.4 mjy and (2) the QSO detection rate grows 
very slowly with improvements in sensitivity. In the 
flux-density range 2.4 < 5(mJy) < 1000, the unbinned 
maximum-likelihood power- law fit ([Crawford et al.l[l970l) 
to the brightness-weighted differential source count of 
these QSOs is S 2 n(S) = (5.0 ±0.4) x S om±0 - 02 Jy sr" 1 ; 
that is, 

n(S) cx S~P , (5) 

where ft = 1.20 ± 0.02 is much lower than the static Eu- 
clidean j3 — 2.5. The Kolmogorov-Smirnov (K-S) test for 
goodness of fit yields the probability P > 0.33 that the 
source count is consistent with the power-law fit. There 
are no gaps or other features in this QSO flux-density 
distribution to suggest more than one source population 
above S — 2.4 mjy, where we expect primarily AGN- 
dominated radio emission. 

Because the flux-density range spanned by these QSOs 
is large, it is also useful to consider the number nio(S) 
of radio sources in each logarithmic flux-density interval 
centered on S. We define nio by 

n 10 (S)d log S = n{S)dS (6) 

so that 

n w (S) = ]n(10)£n(S) oc S 1 ' . (7) 

For our 0.2 < z < 0.45 QSOs, (1-/3) = -0.20±0.02. The 
QSO counts are "flat" in the sense that nio(S) is nearly 
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Figure 2. The peaks in the bri ghtness- weighted count 
S 2 n(S) of all ex tragalactic sources {Condon & Mitchell 119841 ; 
IMitchell fc Uondod 119551) near log[S(Jy)] = -1 and log[S(Jy)] = 
—4 match the expected contributions from AGNs (primarily in 
radio galaxies) and star- forming galaxies, respectively, at typical 
redshifts z ~ 1. The low-redshift SDSS QSO count (filled points 
with error bars) is well fit by the power-law (solid line) whose 
slope also matches the count (open circles) of nearby (z < 0.05) 
UGC galaxies whose radio sources ar e powered primarily by AGNs 
(Condon, Cotton, & Brodcrick 2002). All points in this figure are 
on the same ordinate scale, so the vertical separation between the 
points for all extragalactic sources near the top and SDSS QSOs 
with 0.2 < 2 < 0.45 near the bottom shows how the fraction of all 
radio sources that are in low-redshift QSOs drops from about 1% 
near S = 1 Jy to about 0.01% near S = 1 mjy. The nearly equal 
densities of SDSS QSOs and UGC galaxies is only a coincidence. 
Abscissa: log 1.4 GHz flux density (Jy). Ordinate: Differential 
source count n(S) multiplied by S 2 (Jy sr — 1 ). 

independent of S. In each successively lower decade of 
flux density, the expected number of radio detections 
rises by only dex(0.2) - la 60%. This is quite different 
from the "steep" counts of radio galaxies or of all sources 
in a static Euclidean universe, where (1 — 0) ~ —1.5 
regardless of the luminosity function, and the number 
of radio detections should increase by about 3000% per 
decade. 

The brightness- weighted count of our 0.2 < z < 0.45 
QSOs clearly peaks at the highest flux densities (Fig- 
ure [2]). The strongest 12 radio sources (those with 
S > 900 mjy) comprise only 1% of the 1313 QSOs but 
contribute more than half of the total radio flux density 
from the whole sample, even if all 1148 radio sources not 
in the NVSS catalog had flux densities equal to their 
5u = 2.4 mjy NVSS upper limits. 

2.3. The l.Jf. GHz Luminosity Function and its 
Evolution, 0.2 < z < 0.45 

The 1.4 GHz spectral luminosity distribution of QSOs 
in our sample is shown in Figure [3] as a function of the 
comoving volume in Q = 2.66 sr between the redshift 
cutoff at z = 0.2 and the redshift of each QSO. The QSO 
sample is strictly volume limited and the radio detec- 
tions are complete above the S = 2.4 mjy NVSS limit 
indicated by the curved line. Consequently the plotted 
density of points above the curved line is proportional to 
the 1.4 GHz spectral luminosity function 

(8) 
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Figure 3. The 1.4 GHz luminosity distribution of QSOs detected 
by the NVSS is shown as a function of the volume enclosed be- 
tween the redshift of each QSO and the sample's lower redshift 
limit z = 0.2. The S = 2.4 mjy NVSS sensitivity limit is indi- 
cated by the curve. The QSO sample is volume limited, so the 
density of points above the curve is proportional to its 1.4 GHz 
luminosity function. Symbols distinguish QSOs by their absolute 
magnitudes: large filled symbols for M\ < —24, small filled symbols 
for -24 < Mi < -23.5, and open symbols for -23.5 < Mi < 23. 
Abscissa: comoving volume V (Gpc 3 ). Ordinate: log 1.4 GHz 
spectral luminosity. Parameter: redshift z. 



where m = dex(0.4) = 1 "magnitude" and p(L u )dL v 
is the comoving space density of sources in the spectral 
luminosity range L v and L u +dL u . The QSOs in Figured] 
have been grouped into three absolute-magnitude bins: 
Mi < -24 (large filled symbols), -24 < Mi < -23.5 
(small filled symbols), and —23.5 < M; < —23 (open 
symbols) . 

Figure[3]shows that: (1) The 1.4 GHz spectral luminos- 
ity function has no statistically significant features (re- 
gions of significantly low or high density) in the sampled 
range 24 < log[I/i.4 gh~z(WHz -1 )] < 27 at any redshift 
0.2 < z < 0.45. (2) The radio and optical luminosities 
do not appear to be correlated. However, the radio lu- 
minosities span three decades, while only five of the 163 
QSOs are brighter than Mj = —25.5, one decade brighter 
than the Mj = —23 optical cutoff, so a linear correlation 
would be difficult to discern even if there were one. The 
fact that a linear correlation would be difficult t o discern 
sugge sts that studying radio/optical ratios R (|SchmidtJ 
119701 ) instead of flux densities S is not likely to yield new 
insights i nto low-redshift opti cally selected QSOs. We 
note that iMahonv et al.l (|2012[ ) also found no correlation 
of radio and optical luminosities in their sample of low- 
luminosity X-ray selected QSOs. (3) The plotted density 
of points increases to the right, indicating radio density 
evolution. 

We calculated the 1.4 GHz spectral luminosity function 
of color-selected QSOs brighter than M; = —23 in the 
redshift range 0.2 < z < 0.45 using the standard 1/V max 
method. If there are N sources in the logarithmic bin of 
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Figure 4. The 1.4 GHz luminosity function of QSOs in the red- 
shift range 0.2 < z < 0.45 based on NVSS detections is shown 
by data points at logfLi^ GHz(W Hz -1 )] = 23.6 to 27.2 by 0.4 
and an upper limit at 27.6. The c onti nuo us c urve approximat- 
ing the QSO data follows equations 1111 and 1131 The dotted line 
is a power-low extrapolation to lower luminosities, cut off below 
logfLx 4 GHzfWHz -1 )] ~ 20 to ensure that the number of radio 
sources equals the number of QSOs. Abscissa: log 1.4 GHz spec- 
tral luminosity (W Hz) -1 . Ordinate: log comoving space density 
(Mpc — 3 ) of sources per "magnitude" = dex(0.4) luminosity range. 



width m centered on spectral luminosity L v , then 



^) = e(f-) 

i=i \ ^ max / 1 



(9) 



where V max is the comoving volume in the sample solid 
angle £1 = 2.66 sr between the minimum redshift z = 0.2 
and the maximum redshift z max < 0.45 at which the 
source could have been detected. The rms statistical 
uncertainty of p m (L„) is 



E 



i 



v m 



1/2 



(10) 



The binned luminosity function of QSOs averaged over 
the redshift range 0.2 < z < 0.45 is shown in Figure 01 
The data points can be fit by the power law 



log[p m (mag 1 Mpc 3 )] 



4.21 



-0.mog[L 1AGliz (WRz- L ) 



(11) 



in the luminosity range 23.4 < log[Li.4GHz (WITT 1 )] < 
25.8. The power-law slope 



d\og(p m )/d\og(L v ) 



-0.16 



(12) 



is close to the slope (1— ff) ~ —0.2 of nio(S), as suggested 
in Section |2~21 It is somewhat flatter than the luminosity 
function of UGC radio galaxies powered by AGNs (upper 
curve in Figure [4j. 

At higher luminosities, the luminosity function must 
fall below this flat power law; otherwise the spectral 
power density contributed by the most luminous sources 
would diverge. The amount Alog(p m ) by which the ob- 
served luminosity function falls below the power-law fit 
can be approximated by the (somewhat arbitrary, but 
simple and yielding a satisfactory fit to the limited data) 



quadratic form indicated by the lower continuous curve 
in Figure |4] 



Alog[p m (mag x Mpc 3 )] = 
[(log[ J L 1 . 4G Hz(WHz- 1 )] - 25.8)/1.6] : 



(13) 



when log[Li. 4G Hz(WHz x )] > 25.8. Finally, the total 
number of optically selected QSOs is finite (N = 1313 
over f2 = 2.66 sr in the redshift range 0.2 < z < 0.45), 
so the total number of radio sources predicted by the 
luminosity function should not exceed this number. Any 
extrapolation of the fitted luminosity function to lower 
luminosities (dotted line in Figure |4j must cut off below 
log[ii.4 ghz(W Hz - )] ~ 20 prevent the total number of 
radio sources from exceeding the number of QSOs. 

The NVSS catalog should contain all sources with spec- 
tral luminosities log[Li.4GHz(W Hz -1 )] > 24.2 through- 
out the redshift range 0.2 < z < 0.45. Low-redshift radio 
sources this strong are probably dominated by emission 
from AGNs, not starbursts. The evolution of their lu- 
minosity function with redshift may be characterized in 
terms of pure density evolution (shifting the luminosity 
function horizontally in Figure 0]) 



\og[p m (L,z)} =lag[g(z)p m (L,0)] 



(14) 



or pure luminosity evolution (shifting the luminosity 
function vertically in Figure |4]) 



\og{p m (L,z)]=\og{p m (L/f(z),0)} , 



(15) 



where the functions g(z) and f(z) specify the amounts 
of density and luminosity evolution, respectively. For 
power-law density evolution g(z) = (1 + z) 5 , 5 = +6.0 ± 
1.8 is needed to match the 1.4 GHz data in this portion of 
the redshift-luminosity plane, and the fit is satisfactory 
[P(x 2 ) > 0.2]. Thus the evolving 1.4 GHz luminosity 
function of QSOs in the redshift range 0.2 < z < 0.45 is 
well approximated by replacing Equation 1111 with 

log^mag^Mpc -3 )] = -5.06 
-0.161og[ J L 1 .4GHz(WHz- 1 )] +6.01og(l + z) (16) 

and leaving the cutoff Equation [13] unchanged. 

Pure luminosity evolution of the form / (z) = (1 + z) x 
is indistinguishable from pure density evolution in the 
power-law portion of a luminosity function. However, 
the small power-law slope e « —0.16 of the QSO luminos- 
ity function means that very strong luminosity evolution 
is needed to have the same effect as moderate density 
evolution: A ~ —S/e ~ 38. The high-luminosity cut- 
off cannot evolve this rapidly because there are only two 
sources with log[Li.4GHz(W Hz -1 )] > 26, so pure lumi- 
nosity evolution does not fit our radio data. 

2.4. Why is the QSO Flux-density Distribution So Flat? 

The flat flux-density distribution n\o(S) of radio 
sources in low-redshift optically selected QSOs does 
not prove that they are significantly different from 
the radio sources in nearby optically selected galaxies. 
The open circles in Figure [5] show the weighted count 
S 2 n(S) of AGN-dominate d radio sources in low-redshift 
iz < 0.05) UGC galaxies (|Condon. Cotton, fc Broderickl 
2002); their power-law slope is almost the same. (Their 
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sky density is also nearly the same, but that is just a co- 
incidence that depends on the size of the UGC sample.) 

Why do these optically selected QSO and galaxy sam- 
ples have such flat 1.4 GHz flux-density distributions? 
Consider the flux-density distributions in two extreme 
cases: (1) In an infinite static Euclidean universe, the 
slope of nio(S) would be 1 — j3 = —1.5, for any luminos- 
ity function. (2) If all of the sources in a sample were 
at exactly the same distance, the slope of nio(S) would 
equal the slope of the luminosity function p m (i). For the 
more realistic case of a sample of sources spanning a finite 
rcdshift range z m - m to z max , the flux-density distribution 
will be the convolution of the luminosity function with a 
smoothing function whose width is roughly (z max /^ m i n ) 2 , 
the range of flux densities for a source of fixed luminosity 
as it is moved across the redshift range. 

For both the SDSS QSO and UGC galaxy samples, the 
power-law portions of their 1.4 GHz luminosity functions 
(Figured]) are much wider than (z ma _ x / z m - m ) 2 ~ 6 for 
steep-spectrum QSOs confined to 0.2 < z < 0.45. About 
7/8 of the volume in < z < 0.05 containing the AGN- 
dominated UGC galaxies is in the narrow redshift range 
0.025 < z < 0.05 for which the smoothing width is only 
w 4. In the small-smoothing limit, the power-law slope 
(1 — P) of mo naturally approaches the power-law slope 
dlog( j oio)/(ilogL l/ of the differential luminosity function 
expressed as a source density per log radio luminosity. 
Thus the flat and nearly equal flux-density distributions 
of both the SDSS QSO and the AGN-dominated UGC 
galaxy samples are telling us only that both samples 
have relatively broad and flat radio luminosity functions. 
If the AGN-powered radio sources in both optically se- 
lected samples are intrinsically similar, we should expect 
our QSOs to have a distinctly sub-Euclidean flux-density 
distribution. 

The observed (3 « 1.2 is also consistent with 
relati vistic beaming, whi c h tri es to make /3 » 
4/3 dScheuer fc Readheadl I1979Q . However, rela- 
tivistic beaming of the modestly luminous {typical 
log[Li.4GHz(WHz) _1 ] ~ 25} NVSS radio sources in our 
low-redshift QSO sample is probably insufficient to in- 
fluence the radio flux-density distribuion. In beaming 
models, the most luminous sources are the most highly 
beamed, so that most of their flux densities arise in highly 
relativistic (7 ~ 5-10) components approaching within 
~ 7 -1 radians of our line of sight. Relativistic jets are 
compact and one-sided, but the most luminous sources 
in our sample are not. Table [1] lists 26 sources with 
(apparent) log[Li.4GHz(W Hz -1 )] > 26, but only seven 
have even half of their flux density unresolved by FIRST; 
the rest are too extended (linear diameter > 5 kpc) and 
centered on their optical QSO positions to be strongly 
boosted. The seven fairly compact sources contribute 
only 34% of the total apparent luminosity of the 26 lu- 
minous sources, so removing them would reduce the total 
luminosity of the sample by at most Alog(L„) = —0.18. 
What we see is what we get, and we can calculate reason- 
ably accurate intrinsic luminosities and luminosity func- 
tions by assuming isotropic radio emission. 



2.5. NVSS Statistical Detections: Numerous Faint 
Sources Powered by Star Formation? 



To detect radio sources weaker than S — 2.4 mJy, we 
originally measured NVSS peak flux densities at the po- 
sitions of all N = 179 color-se lected QSOs in the n arrow 
redshift band 0.2 < z < 0.3 (jKimball et al.l 1201 lj ). The 
NVSS image noise has a nearly Gaussian amplitude dis- 
tribution with rms a = 0.45 ± 0.01 mJy beam -1 and 
zero mean (jCondon et al.llT99 8). A starburst with lumi- 
nosity log(i/i©) pa 11 obeying the FIR/radio correla- 
tion should produce a radio source with spectral index 
a w —0.7 and spectral luminosity log[L„(WHz _1 )] w 
22.5 at v = 1.4 GHz in the source frame. Such a source 
would yield an observed 1.4 GHz flux density ranging 
from S = 0.12 mJy at z = 0.3 to S = 0.30 mJy at 
z = 0.2. The NVSS is not sensitive enought to detect 
such faint sources individually, but the distribution of 
peak flux densities on QSO positions can be used to de- 
t ect a large (N > 100 ) sample of them statistically. 

Whi te et al.l (j2007l ) discovered a potential problem 
with statistical detections: the nonlinear CLEAN bias as- 
sociated with the high dirty-beam sidelobes of snapshot 
surveys like FIRST and NVSS can produce a "snapshot 
bias" that lowers the image peak flux densities of sources 
fainter than the catalog limit. The snapshot bias for a 
faint source is about 40% of the peak flux density on a 
FIRST image. To see if snapshot bias might affect the 
peak flux densities of faint QSOs on NVSS images, we 
followed their technique and measured NVSS peak flux 
densities on the positions of faint sources found in the 
sensitive (a « 0.023 mJy beam -1 ) 1.4 G Hz VLA image 
cover ing the Spitzer First-Look Survey (jCondon et al.l 
2003). This image was made from data having very 
good (u, u)-plane coverage, so its snapshot bias should 
be extremely low. The results of this comparison are 
listed in Table O They show that the NVSS snapshot 
bias is too small to measure over the flux-density range 
0.115 < 5p < 0.460 mJybeam -1 . 

Figure shows the distribution of NVSS peak flux 
densities S p at the optical positions of the 179 QSOs 
with 0.2 < z < 0.3 (continuous histogram) and at the 
179 "blank sky" positions exactly one degree north of 
the QSOs (dotted histogram). The blank-sky distribu- 
tion is well fit by a Gaussian with mean (S p ) = 0.00 ± 
0.04 mJy beam" 1 and rms a = 0.48 ±0.04 mJy beam -1 
(dotted curve) consistent with NVSS image fluctuations. 
The whole distribution of peak flux densities on the QSO 
positions is clearly displaced from zero and has a long 
positive tail, indicating that the NVSS has detected most 
0.2 < z < 0.3 QSOs statistically at typical peak flux den- 
sities S p ~ 0.3 mJy beam -1 , which corresponds to an av- 
erage spectral luminosity log[Li.4GHz(W Hz -1 )] ~ 22.7 
in this redshift range. 

Our high detection rate is not consistent with a power- 
law extrapolation of the flat flux-density distribution of 
AGN-dominated QSOs, which would yield the dashed 
curve in Figure [2] The high detection rate of faint 
sources requires a peak or "bump" in the luminosity func- 
tion near log[Zi.4GHz(W Hz -1 )] ~ 22.7. This bump is 
confirmed by individual detec tions of nearly all o f the 
0.2 < z < 0.3 QSOs at 6 GHz (jKimball et alJlMl . 

We hypothesize that the bump in the QSO luminos- 
ity function corresponds to the emergence of the host 
starburst galaxies as the dominant contributors to the 
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Figure 5. The continuous histogram shows the distribution of 
NVSS peak flux densities for the volume-limited sample 179 SDSS 
QSOs in the redshift range 0.2 < z < 0.3. The dotted histogram 
(colored blue in the online version) indicates the distribution of 
NVSS peak flux densities at "blank sky" positions one degree north 
of the QSOs, and it is well fit by a Gaussian of mean (S p ) = 
-0.01 ± 0.04 mjy beam -1 and rms a = 0.48 i 0.04 mjy beam 
(dotted curve, colored blue in the online version). The dashed curve 
(colored green in the online version) is the distribution predicted by 
extrapolating the observed luminosity function of SDSS QSOs to 
lower luminosities, and the continuous black curve represents the 
distribution predicted by our best-fit model in which both AGNs 
and star-forming host galaxies contribute to the radio emission of 
QSOs. 

radio emission. Figure [6] compares the 1.4 GHz lumi- 
nosity functions of low-redshift , (z < 0.05) UGC galaxies 
([Condon. Cotton, fc Broderi ck 2002{) whose radio emis- 
sion is powered primarily by AGNs (continuous green 
curve) or by recent star formation (continuos red curve) 
with a simple model for the radio luminosity function 
of 0.2 < z < 0.3 QSOs. The continuous black curves 
in Figures [5] and |5] correspond to our model luminosity 
function that fits both the NVSS detections and the S p 
distribution of nondctcctions. Our model assumes that 
the total radio emission of each QSO is the sum of sta- 
tistically independent AGN- and starburst-powered com- 
ponents. The lower green curve in Figure [5] is similar to 
the extrapolated luminosity function of Figure 21 con- 
strained by NVSS detections above log[L(WHz -1 )] w 
24, constrained statistically in the luminosity range 
22 < log^WHz- 1 )] < 24 (dashed green line), and at 
lower luminosities by requiring that the number of ra- 
dio sources equal the number of QSOs (dotted green 
curve). The dashed red curve is the radio luminosity 
function of sources powered by recent star formation in 
the host galaxies needed to match the observed total lu- 
minosity fu nction. The direct de tections of faint sources 
reported in iKimball et al.l ()2011[ ) match the dashed red 
curve very well. It suggests that most QSOs host star- 
bursts whose luminosities are comparable with the more 
luminous starbursts in nearby galaxies. The typical 
star-formation rate can be estimated from the peak of 
the dashed red curve: log[Li.4 ghz(WHz)] ~ 22.7 cor- 
responds to a star-formation rate SFR ~ 20M Q yr _1 
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Figure 6. Comparison of the 1.4 GHz luminosity functions of 
nearby galaxies whose radio sources are powered primarily by star 
formation (continous red curve) or by AGNs (continuous green 
curve) with our model for the luminosity function of0.2<z<0.3 
QSOs (black curve). The data points are from Figure [4] with 
log[p m (Mpc ~ 3 m ag -1 ) lowered by 0.20 to correct for density evo- 
lution (Sec. 12.311 . The green dashed line is the power-law extrap- 
olation of the AGN-only contribution, and the green dotted curve 
shows how this extrapolation might be truncated so the number of 
AGNs does not exceed the total number of QSOs. The red dashed 
curve shows the luminosity function of the starbursts only. The 
black curve is the convolution of the AGN and starburst curves. 

(|Condonlll992T ). 

2.6. Stacking 

We note that measuring the full flux-density distribu- 
tion of undetected QSOs is not the same as "stacking," 
which collapses the distribution to a single statistic, the 
mean value of the peak flux density (S p ) of the TV QSOs 
fainter than the catalog limit Su m . The implicit assump- 
tions needed to justify the use of stacking are (1) (S p ) is 
a good statistic for describing the QSO flux-density dis- 
tribution and (2) stacking divides the error in this mean 
by about N 1 / 2 . These assumptions fail if the flux-density 
distribution is a nearly flat power law, as we found for 
the cataloged NVSS QSOs. When nio(S) is nearly in- 
dependent of S, then n(S) ~ S -1 and (S'p) is domi- 
nated by the small number 77 <C N of sources just fainter 
than Sum, the fractional statistical uncertainty in (S'p) is 
about (Sp)/?] 1 / 2 > (SpVV 1 / 2 , and (S p ) depends more 
sensitively on the survey limit Sn m than on the intrinsic 
properties of the optically selected QSO sample. 

Consider a QSO sample whose differential flux-density 
distribution is close to n(S) oc S _1 over a wide (S m - ln <C 
S max ) flux-density range S m - m to S max . The value of S min 
is too low to measure; it is only inferred from the fraction 
/ of detected QSOs and the integral constraint that all 
QSOs must be stronger than S min . The average flux 
density (S) of unresolved sources (so S p ~ S) in a stack 
of images with S < Sn m is 



(S) 



Su 



S n 



fs*2 Sn(S)dS 
/J 1 "" n(S)dS ~ ln(S lim /S min ) 



(17) 
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If the fraction / of individually detected sources is not 
close to unity, then Sn m 3> S m ; n and the numerator 
Slim — Smin ~ ftim- The number of sources per log- 
arithmic flux-density interval is constant between S m i n 
and S max when n(S) oc S -1 so 



In 



Sri 



Sli. 



/In 



Using 



yields 



In | ^1 =ln 



In 



(S P ) 



Mm 



Si; 



lim 



Sr, 



(l-/)ln(Smax/S n 



(18) 



(19) 



(20) 



The stacking flux density (5 P ) depends almost linearly 
on the detection limit Sn m and only logarithmically on 
the intrinsic QSO flux-density distribution. For exam- 
ple, the NVSS detected / = 37/179 « 0.21 of the DR7 
QSOs with 0.2 < z < 0.3 above S lim = 2.4 mJy If 
Smax/Smin ~ 10 8 , then Equation l20l predicts (S p ) » 
Sum/ 15 ~ 0.16 mJy beam -1 , half of which is contributed 
by the small expected number n ks 5 of sources with 
Siim/2 < S < gw 



Whi te et al.l ( 20071 ) made numerical simulations of 
stacking with Gaussian and exponential distributions 
and also concluded that the stack mean does not pro- 
vide a very robust measurement. They studied the stack 
median as well, which should be better than the mean 
because it is less sensitive to the strongest sources. How- 
ever, their simulations showed that the stack median re- 
covered from exploiting low-SNR data (the usual reason 
for stacking) approaches the true mean, so they inter- 
p reted their stack med ians as means. 

iMahonv et al.l (|2012f l has independently noted several 
limitations of stacking. 

We believe it is both safer and more sensitive to fit 
a model to the entire amplitude distribution (Figure [5]) 
than to rely on any stacking statistic. In particular, if the 
whole S p distribution is shifted to positive flux densities, 
we can conclude that most of the QSOs are detectable 
radio sources. 

2.7. Evolution of the Faint Radio Sources, 
0.2 < z < 0.45 

Encouraged by our statistical detection of 1.4 GHz 
emission at starburst luminosities from most of the 179 
QSOs with 0.2 < z < 0.3 and its con firmation by individ- 
ual detections ([Kimball et al.l[20TTI ). we measured NVSS 
peak flux densities at the positions of all 1313 QSOs with 
0.2 < z < 0.45 and divided them into five redshift bins 
(0.2 < z < 0.25, 0.25 < z < 0.3, 0.3 < z < 0.35, 
0.35 < z < 0.4, and 0.45 < z < 0.5) to constrain the 
evolution of the hypothesized starburst radio emission 
over lookback times from 2.4 to 4.6 Gyr. The results 
are shown as histograms in Figure [7] The histogram 
of peak flux densities in 1313 "blank sky" positions one 
degree north of the QSOs is plotted in the top panel, 
and the matching smooth curve shows the Gaussian fit 
with mean (S p ) = — 0.004 ± 0.013 mJy beam -1 and rms 
a = 0.458 ± 0.013 mJy beam -1 . The lower five panels 




•10 12 
1.4 GHz S (mJy/beam) 



Figure 7. Histograms of NVSS peak flux densities for 1313 
SDSS QSOs in five redshift ranges spanning 0.2 < z < 0.45 fit- 
ted with an evolutionary model (smooth curves). The top panel 
shows the distribution of peak flux densities on "blank sky" refer- 
ence positions one degree north of the QSOs; it is nearly Gaus- 
sian with mean (S p ) = —0.004 ± 0.013 mJy beam -1 and rms 
a = 0.458 ± 0.013 mJy beam -1 . 

contain the histograms of NVSS peak flux densities at 
the QSO positions, separated into redshift bins. 

The NVSS detection rate in the 0.2 < z < 0.45 QSO 
sample is only 163/1313 « 12% above 2.4 mJy, so 88% of 
the 1313 QSOs have Sp < 2.4 mJy beam -1 and models 
that match the numbers of sources in each redshift bin 
are constrained primarily by the numbers and hence den- 
sity evolution of the optically selected QSOs. The radio 
data determine the flux-density distribution within each 
redshift bin and hence constrain the luminosity evolution 
of these low-luminosity radio sources. 

As in Section 12.31 we describe the density evolution of 
the optically selected QSOs by p m (z) oc (1 + z) s . The 
best fit, shown by the smooth curves in Figure [71 is for 
S = +8.5 ±2. The best fit is poor [P{x 2 ) - 0.01] because 
the numbers of optically selected QSOs in the redshift 
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bins 0.2 < z < 0.25 (N = 85) and 0.25 < z < 0.3 
(N = 94) are nearly equal but the comoving volumes 
(0.402 Gpc 3 and 0.572 Gpc 3 , respectively) are not. These 
numbers favor slightly negative density evolution (<5 < 0) 
in the 0.2 < z < 0.3 redshift range. The positive value 
5 = +8.5 ± 2 for the faint sources over the whole range 
0.2 < z < 0.45 is consistent with the value S = +6.0±1.8 
for the stronger (primarily AGN-powered) radio sources 
with 0.2 < z < 0.45. 

If we model 1.4 GHz luminosity evolution of the fainter 
radio sources by the form L v (z) oc (1 + z) x , the best fit 
is consistent with no luminosity evolution: A = — 1 ± 3. 

All of these results are actually consistent with each 
other and can be summarized as follows: The optically 
selected QSOs undergo strong density evolution over the 
reshift range 0.2 < z < 0.45. Both the weak and strong 
radio sources undergo a comparably strong density evo- 
lution with no detectable luminosity evolution. The like- 
lihood of radio emission at any spectral luminosity from 
any particular QSO is independent of redshift in the 
range 0.2 < z < 0.45. Although the space density of op- 
tically selected QSOs evolves dramatically, the likelihood 
and strength of radio emission from each QSO appears 
to be independent of redshift. 

3. THE MAGNITUDE-LIMITED 1.8 < Z < 2.5 QSO SAMPLE 

Many historical QSO samples favor the redshift range 
1.8 < z < 2.5 in which the Lya line is redshifted to 
the blue op tical band. To te st the statistical signifi- 
cance of the iMiller et al.l ([1990) bimodal radio luminos- 
ity distribution of such high-redshift QSOs, we chose a 
magnitude-limited sample of all 2471 color-selected DR7 
QSOs brighter than m T = 18.5 in that redshift range. 
These QSOs are all extremely luminous because our red 
magnitude limit corresponds to M\ « —26.9 at z = 1.8 
and Mi w -27.7 at z = 2.5. 

The NVSS detected radio emission stronger than S = 
2.4 mJy from only 191 (8%) of them (Table [3]), consis- 
tent with the usual low detection rates at this sensitiv- 
ity. The brightness- weighted 1.4 GHz flux-density dis- 
tribution S 2 n(S) of the NVSS detections is shown by 
the points with error bars in the lower part of Figure [H 
The power-law fit to their brightness-weighted differen- 
tial source count is 



S 2 n(S) = (12.9 ±1.0) x S 



0.99±0.02 



Jy sr 



(21) 



in the flux-density range 2.4 < 5(mJy) < 1000. Once 
again, the QSO counts are extremely flat in the sense that 
the number of sources per decade of flux density is nearly 
constant: Sn(S) cx 7110(5*) oc S~ 001 . The strongest 13 of 
the 2471 QSOs (« 0.5%) account for more than half of 
the total radio flux density from the high-redshift QSO 
sample. 

Figure |H] shows the luminosity distribution of these 
QSOs as a function of the comoving volume beyond z = 
1.8 in n = 2.66 sr. Our magnitude-limited high-redshift 
QSO sample is volume limited for those QSOs brighter 
than Mi = -27.7 (filled points) above the S = 2.4 mJy 
line, where the plotted density of filled points is pro- 
portional to the evolving radio luminosity function. In 
the redshift range 1.8 < z < 2.5 the density evolution 
of radio sources stronger than log[L(WHz -1 )] = 26 in 
high-redshift QSOs with Mi < -27.7 is small or slightly 
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Figure 8. The peaks in the brightness-weighted count of all extra- 
galactic sources near log[S(Jy)] = —1 and log[S'(Jy)] = —4 match 
the expected contributions from AGNs (primarily radio galaxies) 
and star- forming galaxies, respectively, at typical redshifts 2 ~ 1. 
The high-redshift SDSS QSO count (filled points with error bars 
in the lower part of the plot) of sources stronger than the NVSS 
detection limit S = 2.4 mJy is well fit by a power-law (solid line) 
for which the slope of S 2 n(S) is +0.99 ± 0.02. The dashed and 
dotted curves illustrate the range of faint-source counts needed to 
match the NVSS peak flux-density distribution on the positions of 
high-redshift QSOs fainter than 2.4 mJy. Abscissa: log 1.4 GHz 
flux density (Jy). Ordinate: Differential source count multiplied 
by S 2 (Jy sr" 1 ). 

negative, as expected: if p m oc (1 + z) s , 8 = —0.3 ± 1.8. 
The correlation of radio and optical luminosities is weak; 
if anything, the most luminous radio sources are under- 
represented among the optically most luminous QSOs 
(large filled points) , many of which have radio luminosi- 
ties near log[L(WHz -1 )] = 26. As expected from the 
1.4 GHz source count, the 1.4 GHz spectral luminos- 
ity function averaged over the 1.8 < z < 2.5 redshift 
range is nearly flat (Figure [TO]) until it cuts off around 
log[L l/ (WHz -1 )] » 28.8, the highest radio luminosity 
found in this huge (77 Gpc 3 ) comoving volume. 

3.1. NVSS Statistical Detections: Faint Sources 
Powered Primarily by Star Formation? 



Using the technique described in Sec. 12.51 we ob- 
tained the NVSS peak flux densities at the positions 
of all 2471 high-redshift QSOs and at the 2471 "blank 
sky" positions offset by 1 deg of declination. Figure [IT] 
shows the observed distributions of peak flux densities 
on the QSO positions (continuous histogram) and on 
the blank sky points (dotted histogram). The contin- 
uous histogram is clearly shifted to the right, indicat- 
ing that most of the QSOs are weak radio emitters, 
and their median peak flux density is (S) « 0.05 ± 
0.01 mJy beam -1 . The blank sky distribution is well 
approximated by a Gaussian (dotted curve) with mean 
(S p ) = -0.014±0.009 mJy beam" 1 and rms a = 0.459± 
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Figure 9. The 1.4 GHz luminosity distribution of QSOs detected 
by the NVSS is shown as a function of the volume enclosed be- 
tween the redshift of each QSO and the sample's lower redshift 
limit 2 = 1.8. The S = 2.4 mjy NVSS sensitivity limit is indicated 
by the curve. Symbols distinguish QSOs by their absolute mag- 
nitudes: large filled symbols for M\ < —28.7, small filled symbols 
for -28.7 < M; < -27.7, and open symbols for -27.7 < M ; . The 
QSO sample is volume limited for QSOs with M\ < -27.7 (filled 
points), so the density of filled points above the curve is propor- 
tional to the 1.4 GHz luminosity function of QSOs more luminous 
than M[ = —27.7. Abscissa: comoving volume V (Gpc 3 ). Or- 
dinate: log 1.4 GHz spectral luminosity. Parameter: redshift z. 
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Figure 10. The 1.4 GHz luminosity function of QSOs more lumi- 
nous than Mi = —27.7 in the redshift range 1.8 < z < 2.5 is shown 
by data points at log[ii. 4 gh z (W Hz -1 )] = 25.8 to 28.6 by 0.4 and 
an upper limit at 29.0. Abscissa: log 1.4 GHz spectral luminos- 
ity (W Hz) -1 . Ordinate: log comoving space density (Mpc - 3 ) of 
sources per "magnitude" = dex(0.4) luminosity range. 
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Figure 11. The continuous histogram shows the distribution of 
NVSS peak flux densities for the magnitude-limited sample of 2471 
SDSS QSOs brighter than m r = 4-18.5 in the redshift range 1.8 < 
z < 2.5. The dotted histogram (colored blue in the online version) 
indicates the distribution of NVSS peak flux densities at "blank 
sky" positions one degree north of the QSOs, and it is well fit 
by a Gaussian of mean (S p ) = —0.014 ± 0.009 mjy beam -1 and 
rms a = 0.459 ±0.009 mjy beam" 1 (dotted curve, colored blue in 
the online version). The dashed curve (colored green in the online 
version) is the distribution predicted by extrapolating the flux- 
density distribution of NVSS sources powered primarily by AGNs 
to lower luminosities, and the continuous black curve represents 
the distribution predicted by our best-fit model extremely luminous 
star-forming host galaxies may contribute significantly to the radio 
emission of high-redshift QSOs. 

0.009 mjy beam -1 . Simply extrapolating the flat power- 
law flux-density distribution (Equation |2~T|) of sources 
so luminous {log[L( W Hz -1 )] > 25.8} that they must 
be powered primarily by AGNs and convolving that 
flux-density distribution with the blank-sky Gaussian ht 
yields the dashed curve in FigureQT] The dashed curve is 
a poor match to the continuous QSO histogram because 
the extrapolation predicts that the peak flux densities of 
most QSOs are extremely low (< 0.001 mjy beam -1 ). 
It is necessary to invoke a rise or peak in the flux-density 
distribution near S ~ 0.05 mjy (R ~ 0.3) to explain the 
high detection rate implied by the positive offset of the 
QSO histogram. 

The dotted line in Figure [5] indicates the weighted 
faint-source count from the power-law extrapolation that 
(1) matches with the observed count at 2.4 mjy, (2) satis- 
fies the integral constraint that the total number of radio 
sources must equal the total number of QSOs, and (3) 
best fits the observed flux-density distribution of faint 
sources; it is 



n(S) = OTIS' -1 ' 8 Jy -i sr -i ,S > S„ 



(22) 



and n(S) = below S m [ n = 0.018 mjy. Of course, there 
is no reason to expect a break in the count slope to oc- 
cur at the NVSS catalog limit, 2.4 mjy. At the midpoint 
of the QSO redshift range, z = 2.15, that break corre- 
sponds to a spectral luminosity [logL(W Hz -1 )] ss 25.8, 
which is higher than the spectral luminosities of even the 
strongest starbursts. 
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Any more realistic extrapolation with a continuous 
count and count slope at 2.4 mJy must have a peak in 
S 2 n(S) at lower flux densities to match the statistical 
data, as shown by an illustrative model that yields the 
dashed curve in Figure [5] Convolving the model counts 
with the NVSS image noise distribution yields a good fit 
to the observed distribution of QSO peak flux densities; it 
is shown by the continuous curve in Figure [TTI) . The me- 
dian flux density of our high-redshift QSOs in this model 
is (S) w 0.05 mJy, so the median 1.4 GHz spectral lumi- 
nosity in the source frame (calculated assuming a median 
spectral index a = -0.7) is log[(L„(W Hz -1 )] « 24.1 
at z = 2.15. This is a factor of four higher than the 
1.4 GHz spectral luminosity of Arp 220, and it might 
come from an AGN or a starburst. If it originates in 
a starburst obeying the FIR/radio correlation, an ultra- 
luminous [(L/Lq) ~ 4 x 10 12 ] but not hyperluminous 
[(L/Lq) > 10 13 ] starburst is required. 

4. IS QSO RADIO EMISSION BIMODAL? 

The first published claim of a bi modal flux-density 
distribution (|Strittmatter et al"1 [l980) features a strong- 
source peak too narrow to be consis t ent w ith relativis- 
tic beaming. The iStrittmatter et all (jl980f ) flux-density 
distribution plot is reproduced in Figure 1121 The nar- 
row peak in the quantity A log N/A log S (where N is 
the number of sources stronger than S and A log S = 0.2 
is the logarithmic width of each flux-density bin) near 
log[5(mJy)] = 3 is an apparently significant factor of 
dex(2) ~ 100 higher than the long tail of fainter sources. 
We have inserted the actual numbers of radio sources 
contributing to each logarithmic flux-density bin into 
Figure [12] The peak is produced by only two sources, as 
N goes from 1 to 3, so A log A = log(3) - log(l) w 0.48 
and A log N/ A log S « 0.48/0.2 » 2.4 as plotted. How- 
ever, the quantity A log N/A\ogS is not the usual log- 
arithmic differential source count or flux-density distri- 
bution log A N/ A log S. The thin curve shows the dis- 
tribution of A log Nf A log S that would result from the 
same total number of radio sources (ten) having a con- 
stant log AA/ A log 5* in the plotted flux-density range. 
Thus there is no significant peak in the nearly constant 
logarithmic differential flux-density distribution of these 
Q SOs, which is cons istent with our Equation |2~T1 

Mill er et al.l (|1990l ) made sensitive 5 GHz observations 
of 105 bright QSOs with 1.8 < z < 2.5 selected by their 
emission lines. They found nine with log[L(W Hz -1 )] > 
26 but none with 25 < log[L(WHz -1 )] < 26. All nine de- 
tections have S > 20 mJy, far above their 3<r ~ 0.6 mJy 
detection limit. This is evidence for bimodal flux-density 
and luminosity distributions, but at what statistical sig- 
nificance? The nine detections have flux densities 21, 28, 
28, 51, 76, 98, 105, 111, and 136 mJy. The likelihood of 
no weaker detections above 0.6 mJy depends on how the 
flux densities are binned. The bin limits should not be 
chosen a posteriori to minimize this probability. In bins 
of width dex(l) starting at the detection limit, there are 
N = detections with 0.6 < S(mJy) < 6, N = 4 de- 
tections with 6 < 5(mJy) < 60, and N — 5 detections 
with 60 < 5(mJy) < 600. If the flux-density distribution 
of 1.8 < z < 2.5 QSOs is flat in logarithmic bins (as in 
Equation [21]), there should be equal numbers in each bin. 
The average number of detections per bin is (N) = 3, so 
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Figure 12 . The differentia l flux-density distributions of QSOs as 
plotted by Strittmattcr ct al. ( 1980) with numbers of QSOs per 
bin added. 

the Poisson probability of no detections in the faintest 
bin is Pp(0,3) = cxp(-3) « 0.05. Thus the lack of de- 
tections fainter than 6 mJy in this small sample may just 
be the result of bad luck, which we did not encounter in 
our much larger high-redshift QSO sample. 

We prefer to study evolving QSO luminosity functions, 
which arc intrinsic properties of QSOs, rather than flux- 
density distributions or luminosity distributions, which 
depend on observational parameters. We found that the 
luminosity functions p m (L) of radio sources so power- 
ful that they must be produced by AGN are smooth 
and nearly flat power laws with no peaks or other fea- 
tures. There are "bumps" in the radio luminosity func- 
tions of QSOs, but they appear only at such low ra- 
dio luminosities that the resulting radio flux densities 
are well below 1 mJy. The "bump" luminosities are 
consistent with significant starburst contributions (re- 
call Figure H|: log[L(W Hz -1 )] - 22.7 for our moderate- 
luminosity QSOs at low (0.2 < z < 0.45) redshifts and 
log[L(W Hz -1 )] ~ 24.1, corresponding to an (L/Lq) ~ 
4 x 10 12 ultraluminous starburst, for the most luminous 
QSOs at higher redshifts, 1.8 < z < 2.5. 

5. THE AGN AND HOST GALAXY COMPONENTS OF QSO 
RADIO EMISSION 

The radio emission from a QSO is the sum of contribu- 
tions from its AGN and from star formation in its host 
galaxy. If the luminosities of the two contributions are 
nearly independent of each other at every redshift, the 
spectral luminosity function of QSOs will be the con- 
volution of the AGN and starburst spectral luminosity 
functions. 

5.1. Low-redshift QSOs 

The AGN-powered portion of the luminosity func- 
tion in the redshift range 0.2 < z < 0.45 (Figure [4| 
was estimated from 163 NVSS detections of luminous 
sources. There are only 37 NVSS detections in the 
narrow redshift range 0.2 < z < 0.3, too few for sta- 
tistical accuracy, so we estimated the 0.2 < z < 0.3 
AGN luminosity function from the evolving AGN lumi- 
nosity function in the redshift range 0.2 < z < 0.45. 
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The luminosity function in the range 0.2 < z < 0.3 
is lower than the 0.2 < z < 0.45 luminosity func- 
tion by Alogp = -0.20 (Sec. [273]). This estimate is 
shown by the data points with error bars in Figure [H 
For comparison, the 1.4 GHz local (z < 0.05) luminos- 
ity fu nctions of galaxies (|Condon. Cotton, fc Broderickl 
|2002j) whose radio emission is dominated either by AGNs 
(continuous green curve) or by star formation (contin- 
uous red curve) are also shown. The radio luminos- 
ity function produced by AGNs in local galaxies is flat- 
ter and extends to higher radio luminosities than the 
luminosity function of star-forming galaxies, so local 
AGN-powered radio sources are more common above 
log[Li.4 ghz(WHz~ 1 )] = 23 and star-forming galaxies 
have higher space densities than AGNs at lower lumi- 
nosities. The luminosity log [£1.4 gh~z(WHz -1 )] ps 22 of 
the galaxy M82 is typical for a nearby radio-selected star- 
forming galaxy. 

Even allowing for starburst luminosity evolution, it is 
likely that the radio emission of low redshift (0.2 < z < 
0.45) QSOs stronger than log[Li. 4 ghz(WHz^)] = 24 is 
dominated by AGNs. To estimate the luminosity func- 
tion that would result if all QSO radio emission were 
powered by only by AGNs, even at lower luminosites, 
we made the smooth power-law extrapolation (dashed 
green line in Figure [6j> of the observed QSO luminosity 
function to lower luminosities. This extrapolation must 
be truncated below log [L 1.4 ghz(WHz -1 )] ~ 20, as sug- 
gested by the dotted green curve in Figure EJ to keep the 
number of radio sources from exceeding the total num- 
ber of low-redshift QSOs (1313). The details of the AGN 
luminosity function below log [£1.4 gh z (W Hz -1 )] w 22 
are not well known but don't matter because any lumi- 
nosity function in this range always predicts that most 
of the low-redshift QSOs powered by entirely by AGN 
would be too radio quiet (< 0.02 mJy) to be statistically 
detectable by the NVSS. The dashed curve in Figure [5] 
shows the AGN-only flux-density distribution predicted 
by this model for the 0.2 < z < 0.3 QSOs. The distribu- 
tion peak flux density (S p ) « 0.05 mJy beam -1 is much 

lower than the observed (S p ) ~ 0.33 mJy beam - . The 
AGN-only distribution is distinguishable from the blank- 
sky distribution (dotted curve) only because its height 
has been depressed by the loss of QSOs in the long tail 
of strong sources above 2.5 mJy beam -1 , which are off 
scale on this plot. The failure of the dashed curve to fit 
the data suggests that most of the low-redshift QSOs are 
not in "red and dead" elliptical galaxies with little or no 
ongoing star formation. 

Any QSO luminosity function that is consistent with 
both the NVSS direct detections of sources stronger than 
2.4 mJy and the distribution of peak flux densities in Fig- 
ure [5] must rise sharply just below log[L„(WHz -1 )] 
'24 and fall fairly quickly at lower luminosities lest it 
imply more than 13 1 3 SP SS QSOs with 0.2 < z < 
0.45. Kimba ll et al.l ([201 ll ) proposed that radio emis- 
sion from the star-forming host galaxies is the cause of 
that rise. The peak of the the QSO host-galaxy radio 
luminosity function (dashed red curve in Figure |6|) is 
at log[L„(W Hz -1 )] ss 22.7. If the faint radio sources 
are dominated by emission from the star-forming host 
galaxies, their radio spectral indices between 1.4 GHz 



and 6 GHz should be a ~ —0.7, which is the observed 
value (jKimball et al.ll2011l ). 

As it must, the radio luminosity function of the low- 
redshift QSO host galaxies lies below the radio luminos- 
ity function of all nearby galaxies powered primarily by 
star formation, but it is not much lower at the high- 
est radio luminosities because many low-redshift QSOs 
seem to have high star-formation rates. The radio lu- 
minosity function of QSO host galaxies must fall rapidly 
below log[Li. 4 ghzW Hz -1 ~ 20 to satisfy the constraint 
that its integral not exceed the total number of SDSS 
QSOs in 0.2 < z < 0.3. In this respect, it differs from 
the radio luminosity function of all nearby star-forming 
galaxies (red curve in Figure [5]) , which rises slowly but 
monotonically even at very low luminosities. This differ- 
ence is not surprising: most of the faintest nearby radio 
source s lie in low- mass galaxies with low s tar- formation 
rates (|Condon. Cotton, fc Broderic^[2002h . while QSOs 
are found only in massive galaxies. The real surprise to 
radio astronomers accustomed to comparing QSOs with 
radio galaxies is that massive elliptical galaxies with low 
star-formation rates are so rare in our low-redshift QSO 
sample. 

5.2. High-redshift QSOs 

The high statistical detection rate of 1.8 < z < 2.5 
QSOs reported in Sec. l3.1l shows that most of these very 
luminous (median (Mi) ~ —27.5) high-redshift QSOs 
contain radio sources with median spectral luminosity 
log[Li. 4 ghz(W Hz -1 )] « 24.1. This high detection rate 
is not consistent with a power-law extrapolation from 
higher radio luminosities of the AGN luminosity func- 
tion; there must also be a "bump" in the QSO radio 
luminosity function at high redshifts. As illustrated in 
Figure [H the bump must lie in a fairly narrow flux- 
density range —5 < log[5(Jy)] < —3 to be consistent 
with the distribution of NVSS peak flux densities (Fig- 
ure Hip . Such a bump is consistent with, but does 
not by itself require, radio emission from ultraluminous 
[(L/Lq) ~4x 10 12 ] starbursts in the host galaxies of the 
most luminous QSOs. Should we expect that most lu- 
minous, high-redshift QSOs contain ultraluminous star- 
bursts? 

The presence of a luminous QSO implies a high rate of 
radiativcly efficient accretion onto a supermassive black 
hole. During "quas ar mode" or "cold mode" accretion 
(jCroton et al.l l2006h . cold gas is rapidly fed via an ac- 
cretion disk to the black hole. The reservoir of cold gas 
that feeds the AGN may also produce a burst of star for- 
mation, as suggested by recent observations starbursts in 
the host galaxies of high-luminosity (L^ > 10 37 W) X- 
ray QSOs with z > 2 (|Rovilos et al.ll2012l ) and of infrared 
emission from cold dust in the host galaxies of the fairly 
luminous (( M) -25) QSOs fou n d in the Herschel AT - 
LAS survey ([Bonfield et al.ll2011h . iBonfield et all ()2011[ ) 
associated the total (8 fi < A < 1000 fi) infrared luminos- 
ity -Lir with the star-formation luminosity and found it 
is related to the QSO luminosity Lqso vi a 

L m xL e QSO (l + z)t , (23) 

where = 0.22 ± 0.0 8 and C = 1.6 ± 1. Using the 
FIR/radio correlation (|Condon| I1992D and the median 
1.4 GHz spectral luminosity, we estimate the median 
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starburst luminosity log^m/i©) ~ 11.2 for our low- 
redshift QSOs at typical redshift z « 0.35. The me- 
dian absolute magnitudes of the low- and high-redshift 
QSOs are M; = -23.38 ± 0.03 and M; = -27.48 ± 
0.01, for a QSO luminosity ratio « 44. Inserting 
these values into Equation [23] indicates that our high- 
redshift [z w 2.15) QSOs should have a median star- 
burst luminosity log(i/L Q ) = 12.2 ± 0.4 and a me- 
dian radio spectral luminosity log[Li. 4 gh z (W Hz -1 )] = 
23.7 ± 0.4, only one sigma below our observed median 
log[Li. 4 ghz(W Hz -1 )] = 24.1 for SDSS QSOs in the red- 
shift range 1.8 < z < 2.5. 

6. SUMMARY 

We investigated the 1.4 GHz radio emission from large 
optically selected samples of SDSS DR7 QSOs in two 
redshift ranges: 0.2 < z < 0.45 and 1.8 < z < 2.5. As 
expected, only about 10% of the QSOs are sufficiently 
radio loud to be detected above the S — 2.4 mJy NVSS 
catalog limit. Such sources are so luminous that the bulk 
of their radio emission must be powered by AGNs. The 
1.4 GHz luminosity functions of AGN-powered sources 
both samples (Figures [4] and [TO]) are flat power laws, so 
flat that the number of sources per decade of flux den- 
sity or luminosity is nearly constant. These luminosity 
functions show no features or signs of bimodality that 
might indicate two or more distinct QSO types among 
the - 10% of QSOs that are in the NVSS catalog. 

Extrapolating the flat luminosity functions to lower lu- 
minosities predicts that most of the undetected QSOs 
should be extremely faint radio sources. However, 
the distributions of peak flux densities at QSO posi- 
tions on the NVSS images reveal that most QSOs are 
moderately luminous radio sources, so there must be 
peaks or "bumps" in their 1.4 GHz luminosity func- 
tions in both redshift ranges. QSOs with 0.2 < 
z < 0.3 have a median peak flux density S'p w 
0.3 mJy beam -1 (Figure [5]) and median spectral lumi- 
nosity log[Li. 4 G Hz(W Hz" 1 )] w 22.7; QSOs in the red- 
shift range 1.8 < z < 2.5 have a median peak flux density 
S'p w 0.05 mJy beam -1 (Figure lllj) and median spec- 
tral luminosity log [L 1.4 ghz(W Hz -1 )] rs 24.1. The re- 
ality of the statistical detections in the narrow redshift 
range 0.2 < z < 0.3 was confirmed by individual detec- 
tions of nearly all of these QSOs by the VLA at 6 GHz 
(jKimball et al.l [201 If ) . and the median spectral index of 
these sources is (a) » —0.7, typical of optically thin 
synchrotron emission from either AGNs or star-forming 
galaxies. 

We suggest that the faint radio sources found in most 
bright optically selected QSOs are primarily powered 
by star formation. The median radio luminosities cor- 
respond to star- formation rates M ~ 20M Q yr -1 and 
M ~ 500Af© yr -1 for the moderately luminous ((Mi) w 
-23.4) QSOs with 0.2 < z < 0.45 and the extremely lu- 
minous ((Mi) « -27.5) QSOs with 1.8 < z < 2.5, respec- 
tively. The total infrared luminosities of such starbursts 
can be estimated from the FIR/radio correlation. They 
are log(£ IR /£ Q ) « 11.2 for the 0.2 < z < 0.45 QSOs and 
log(i m /i ) « 12.6 for the 1.8 < z < 2.5 QSOs. Such 
powerful starbursts may be fueled by the same cold gas 
reservoirs that flow into the central supermassive black 
holes that power the QSOs themselves. Our starburst 



interpretation of the low-luminosity QSO radio emission 
is supported by (1) far-infrared detection s of cold dust in 
QSO host galaxies (Bonfiel d et al.|[20TlT ). and (2) agree- 
ment of our radio results with their empirical scaling of 
LrR with QSO luminosity and redshift. 

With cm- wavelength radio data alone, it is difficult 
to test our hypothesis that star formation in the host 
galaxies of most QSOs powers the faint radio sources 
that make the "bumps" in our otherwise flat radio lu- 
minosity functions and source counts. One way would 
be to show that the faint radio sources in QSOs fol- 
low the tight FIR/radio correlation obeyed by nearly all 
starburst galaxies. However, it is necessary to distin- 
guish the FIR emission from a large mass of cool dust 
heated by stars from the mid-infrared (MIR) emission 
from a small mass of warmer dust heated by an AGN 
(M oric et afll201Q[ ). so it is preferable that the the dust 
flux density be measured in the Rayleigh- Jeans tail at 
submillimeter wavelengths, by ALMA for example. Al- 
ternatively, ALMA submillimeter spectroscopy of molec- 
ular and ionic lines sensitive to X-ra y dominated re- 
gions might reveal deeply buried AGNs (jRangwala et al.l 

Eml). 
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Table 1 

NVSS-detected QSOs, 0.2 < z < 0.45 









Redshift 




S 


logLi.4 


RA (J2000) DEC 




z 


Mi 


(mjy) 


(W Hz" 1 ) 


07 54 03.61 


+48 14 28 


.1 


0.2755 


-23.478 


7.3 


24.19 


08 06 44.43 


+48 41 49 


.2 


0.3701 


-23.769 


901.7 


26.57 


08 20 15.61 


+59 42 28 


.4 


0.3676 


-25.461 


5.6 


24.36 


08 33 53.88 


+42 24 01 


.9 


0.2491 


-23.903 


248.6 


25.62 


08 36 58.91 


+44 26 02 


.3 


0.2544 


-24.999 


6.6 


24.07 


08 43 10.79 


+39 53 45 


.1 


0.4036 


-23.765 


2.6 


24.12 


08 43 47.85 


+20 37 52 


.5 


0.2273 


-23.034 


59.7 


24.92 


08 49 40.01 


+09 49 21 


.1 


0.3656 


-23.425 


587.2 


26.37 


08 50 39.96 


+54 37 53 


.3 


0.3673 


-23.076 


166.5 


25.83 


08 55 16.21 


+56 16 56 


.8 


0.4422 


-23.759 


5.3 


24.52 


08 56 32.99 


+59 57 46 


.9 


0.2830 


-23.858 


239.4 


25.73 



Note. — Table \T\ is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 



Table 2 

NVSS "snapshot bias" 



FLS flux density 


Number of 


NVSS bias 


(MJy) 


sources 


(^tjy beam -1 ) 


115-160 


937 


-17+15 


160-230 


821 


+7+16 


230-460 


771 


+29 ± 17 


115-460 


2550 


+3 + 9 



Table 3 

NVSS-detected QSOs, 1.8 < z < 2.5 







Redshift 




5 


logLi.4 


RA (J2000) DEC 


z 


Mi 


(mJy) 


(W Hz" 1 ) 


08 06 20.47 


+50 41 24.4 


2.4565 


-28.38 


16.9 


26.75 


08 11 41.98 


+51 57 10.9 


2.1763 


-27.65 


5.1 


26.12 


08 21 53.82 


+50 31 20.4 


2.1326 


-27.93 


58.5 


27.16 


08 27 54.25 


+33 36 04.2 


1.8056 


-27.04 


31.5 


26.72 


08 33 50.61 


+38 39 22.8 


2.0143 


-27.53 


2.7 


25.76 


08 34 00.05 


+43 01 48.1 


2.3916 


-27.39 


19.0 


26.78 


08 37 22.41 


+58 25 01.8 


2.1010 


-27.89 


690.4 


28.21 


08 45 06.24 


+42 57 18.4 


2.0945 


-27.38 


224.0 


27.72 


08 45 47.19 


+13 28 58.1 


1.8834 


-28.40 


413.4 


27.88 


08 50 51.80 


+15 22 15.0 


2.0183 


-27.85 


520.2 


28.05 



Note. — Table [3] is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 



